Dynamical Screening Effect on Local Two-Particle Vertex Functions 



m 
o 



(N 



-f— > 

CO 



i 

c 

o 
o 



> 

00 

oo 
(N 

o 



Li Huang 1,2 and Yilin Wang 3 

1 Science and Technology on Surface Physics and Chemistry Laboratory, 
P.O. Box 718-35, Mianyang 621907, Sichuan, China 
2 Department of Physics, University of Fribourg, Fribourg CH-1700, Switzerland 

3 Beijing National Laboratory for Condensed Matter Physics, 
and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 

(Dated: March 13, 2013) 

In principle, the electronic Coulomb interaction among the correlated orbitals is frequency- 
dependent. Though it is generally believed that the dynamically screened interaction may play 
a crucial role in understanding the subtle electronic structures of strongly correlated materials, 
hitherto we know very little about it. In the Letter, we demonstrate that within the framework of 
single-site dynamical mean-field theory the local two-particle Green's functions \ an d vertex func- 
tions L are strongly modified by the dynamically screened interaction. Since both x and L represent 
the main ingredients to compute momentum-resolved response functions and to treat non-local spa- 
tial correlations by means of diagrammatic extensions of dynamical mean-field theory, it is urgent to 
reexamine previous results by taking the dynamical screening effect into account. The modifications 
should be very considerable. 

PACS numbers: 71.27.+a, 71.10.-w, 71.15.-m, 71.30.+h 



Over the last years, tremendous progresses have been 
made in the studies of strongly correlated materials by 
using the density functional theory (DFT) combined with 
the dynamical mean- field theory (DMFT). [2,0 The un- 
derlying physics in these materials can be generally cap- 
tured by effective Hamiltonian models which are derived 
in the DFT calculations. Subsequently, these models are 
solved numerically or analytically in the DMFT part to 
extract various one-particle and two-particle quantities, 
which are used to compare with the corresponding ex- 
perimental results if available. 

In classical DFT + DMFT calculations, the electron- 
electron interactions among the correlated orbitals are 
usually modeled by Coulomb interaction matrix which 
is parameterized with Hubbard U. It is hardly surpris- 
ing that the calculated results rely heavily on the style 
of Hubbard £7.0 Usually U is assumed to be static. 
However, the charge density fluctuation in correlated or- 
bitals will induce electric fields, which will be screened 
by the higher-lying or lower-lying states in the systems. 
The screening effect will finally result in a frequency- 
dependent renormalization for the Coulomb interaction. 
Thus, in principle the Coulomb interaction is dynami- 
cal, i.e, U = Since the dynamical screening 
effect may exist in strongly correlated materials widely, 
how to treat it is one of the major challenges in study- 
ing the intriguing properties of these materials. Due to 
many technical difficulties, it has been completely ne- 
glected or empirically taken into account by adjusting 
the static value Uq = U(lu = 0)0 or reducing the effec- 
tive bandwidth of correlated orbitals in most calcula- 
tions. Very recently, numerous efforts have been made 
to explore the dynamical screening effect in strongly cor- 
related materials, but very little is known on the 
consequences of frequency-dependent Coulomb interac- 



tion. So far the only thing we can confirm is that the 
dynamical screening effect may play an essential role in 
understanding the fascinating properties of strongly cor- 
related materials, llj, [l2[ such as the plasmon satellites 



and spectral weight transfer in the photoemission spec- 
tra of SrV03 HI, 14 1 and hole-doped BaFe2 As2 , [IH which 
are less emphasized before. 

The local two-particle vertex functions represent a 
crucial ingredient for the calculations of (momentum- 
dependent) dynamical susceptibilities, such as spin den- 
sity wave, optical conductivity, thermopower, and ther- 
mal conductivity, etc. [l6l - l22j Remarkably, it is argued 
that the local vertex functions can be used to calcu- 
late the electric polarization in correlated insulators [23| 
and to identify the fingerprints of Mott-Hubbard metal- 
insulator transition. [24{ Beyond that, the local two- 
particle vertex functions are also the prerequisites for 
the diagrammatic extensions of the DMFT, such as the 
dynamical vertex approximation (DTA) |26| and dual 
fermion (DF) j2?| approach, which aiming at the inclusion 
of non-local spatial correlations. Therefore, to obtain 
reliable and accurate local two-particle vertex functions 
is one of the major tasks in DMFT calculations all the 
time. However, though many numerical tricks and ap- 
proximations were developed in recent years, 1^-21| it is 
still higher difficulty and heavier workload of performing 
calculations at the two-particle level than the one-particle 
level. As a result it is not astonished that the analysis of 
two-particle quantities is usually restricted to the easiest 
cases of optical conductivity, for which one can neglect 
the vertex correction and retain the bubble term only, 
and then the calculation can be simplified to one-particle 
level eventually. [2-0 

In this Letter, we calculate and analyze the local 
two-particle Green's functions \ and vertex functions L 
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is often used to describe the dynamically screened be- 
haviour for Hubbard U : 



FIG. 1. (Color online) The typical Ohmic models used in the 
present calculations for frequency-dependent Coulomb inter- 
action KeU(uj). The Uo and U BCT are fixed to 8.0 eV and 3.0 
eV, respectively. 



for single-band Hubbard model and strongly correlated 
metal SrVOa by using the DMFT and DFT + DMFT 
approaches [l|-|3| respectively. The dynamical screening 
effects for both the Hubbard model and realistic materi- 
als are considered on the same footing in our calculations. 
We found that the calculated results for both \ and T 
are significantly affected by the dynamically screened U. 
Therefore, it is suggested that the previous results and 
conclusions, concerning with the momentum-resolved re- 
sponse functions and diagrammatic extensions of DMFT, 
etc., should be reexamined and reconsidered carefully, 
provided that the dynamical screening effect is sizable 
and can not be ignored completely. 

We now illustrate the analytical Ohmic model0] which 

I 



ReLT(w) = U 



LO c + UJ 



0J r 



(1) 



where a and uj c are the screening parameters, Uq is the 
static interaction, and the screened interaction can be 
defined as U SCI = Uq — 2au c . Note that a characteris- 
tic feature of U (w) of paramagnetic Ni can be captured 
by this Ohmic model, [j, Q Several typical Ohmic mod- 
els used in the present calculations are shown in FigfT] 
In order to include the dynamical screening effect into 
the DMFT and DFT + DMFT frameworks, [ljl(J the 
hardest obstacle has been the lack of a reliable and effi- 
cient quantum impurity solver for the general impurity 
model with a frequency-dependent Coulomb interaction 
U(uj), Fortunately, it seems that this obstacle has been 
overcome by the recent development of hybridization ex- 
pansion continuous-time quantum Monte Carlo impu- 
rity solver (dubbed CT-HYB) proposed by Werner et 
aL,0, IH where a multi-plasmon Lang-Firsov transforma- 
tion is treated exactly in the context of a hybridization 
expansion algorithm for the general impurity model. 

According to the literatures, [H El the general two- 
particle Green's function \ can be defined as follows: 

Xabcd(T a ,T b ,T c ,T d ) = {TrCa{Ta)c\{T b )Cc{T c )c\{Td)) ■ (2) 

Generally these functions are evaluated in frequency 
space, as a function of two fermionic frequencies [v and 
v') and one bosonic frequency (w). Then the local two- 
particle vertex functions r^, w can be derived directly by 
Xaa /U> through the following equation: 



-$G a {y)G a {v + u)S vv ,S aa , - G a (y)G a (v + u})V£, u G a ,(i/)G ir , (1/ + w) 



(3) 



where j3 is inverse temperature. We calculate the Xaa> u 
and r^> with single-site DMFT and DFT + DMFT 
methods by using the powerful CT-HYB quantum impu- 
rity solver, il, 29| An improved estimator for xVa SU S" 
gested by Hafermann et al. [H, [2l| is adopted to suppress 
the numerical noises and yield accurate data. 

For the sake of simplicity, above all we consider a half- 
filled one-band Hubbard model on the Bethe lattice with 
bandwidth it at /3 = 50.0. The static and screened Hub- 
bard U are 8.0 eV and 3.0 eV, respectively. Then we tune 
the a and w c parameters for the Ohmic model to monitor 
the variational trends for selected one-particle quantities 
and local two-particle vertex functions. 

(a) Dynamical screening effect on spin-spin correlation 
function and local magnetic moment in one-band Hub- 
bard model. Previous studies show that the dynamically 



screened interaction can shift the metal-insulator phase 
boundary, Q change the spectral function near the Mott- 
Hubbard gap edge, and renormalize the quasiparticle 
band structure. [ 1 3ll 1 4| Inspired by these works, we study 
the spin-spin correlation function xss( T ) = (S z (0)S z (t)} 
in the Hubbard model with frequency-dependent U to 
gain more insight into the dynamical screening effect. 
The obtained xss( r )j effective local magnetic moment 
Me = V^XSSi and quasiparticle weight Z are shown in 
Figl^l As is seen in this figure, the xss{i~), M e , and Z 
are strongly modulated by the screening frequency pa- 
rameter lo c of the Ohmic model. Obviously, there ex- 
ists a critical oj c separating the metal and Mott insulator 
phases. We just focus on the following two cases: (i) 
The screening frequency lj c is small {uj c < 4.0). The 
Xss(t) approaches a non-zero constant at large times 
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FIG. 2. (Color online) Magnetic moment collapse and metal- 
insulator transition in half-filled one-band Hubbard model 
induced by the dynamical screening effect, (top) Spin-spin 
correlation function xss(t) = (S z (0)S z (t)) . (bottom) Quasi- 
particle weight Z and local spin magnetic moment x/T^ss as 
a function of uj c . The Mott metal- insulator transition zone is 
highlighted by a pink vertical bar. 



which characterizes a frozen moment phase. M e only de- 
creases slightly from its initial value (0.98 \ib at cj c = 0.0) 
under the increment of uj c . And the quasiparticle weight 
Z approximates zero, which means the system should ex- 
hibit completely insulating properties, (ii) The screening 
frequency uj c is large (oj c > 4.0). The behavior of xss(t) 
manifests that the system still retains the frozen moment 
phase. Whereas, M e decreases quickly and monotonously 
from 0.96 to 0.43 i.e, the effective local magnetic 
moment appears apparent collapse. On the other hand, 
the quasiparticle weight Z becomes more and more con- 
siderable when u> c increases, which means that the sys- 
tem goes into a fully metallic state. Further analysis 
on the imaginary part of the calculated low-frequency 
Matsubara self-energy function ImS(io;) indicates that 
it exhibits noticeably fractional power-law behavior with 
respect to iu, which is the signature of the non- Fermi- 
liquid metallic phase. In Fig(2l we observe a sudden drop 
for M e and rise for Z around uj c = 4.0, which corresponds 
to a typical first-order transition. 



FIG. 3. (Color online) Vertex functions vs. the two fermionic 
frequencies v = (2n + l)n//3 and v' — (2n + 1)tt/P as a func- 
tion of the screening frequency cj c for half-filled single-band 
Hubbard model, (top) Real part of spin-up-up component for 

the local two-particle Green's functions: Rex^f w . (bottom) 
Real part of spin-up-up component for the local two-particle 
vertex functions: ReF^ Note that the bosonic frequency uj 
is fixed to be zero. Instead of the absolute values of the Mat- 
subara frequencies just the corresponding indexes are given 
so as to improve the readability of these figures. 



(b) Dynamical screening effect on local two-particle 
vertex functions in one-band Hubbard model. We then 
concentrate our attentions to the local two-particle 
Green's functions Xaa^ an( i ver tex functions 
In Fig|3J only the contour maps for their spin-up-up com- 
ponents (the bosonic frequency u> = 0) are displayed. As 
is clear in this figure, if the dynamical screening effect 
is included in the DMFT calculations, not only the one- 
particle quantities but also the two-particle quantities are 
modified visibly, which are not taken seriously before. 
Let's look at the Rex++ u part (see top panel in FigJ3]) 
at first. The intensity plots of it can be divided into 
two kinds, (i) For uj c < 4.0, one observes a gradual en- 
hancement at the low-frequency zone (central area). But 
neither the singularity nor obvious geometric structures 
is detected, (ii) For uj c > 4.0, a numerical divergence 
emerges at v — > and v' — > zone, where the intense 
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white-black color coding indicates alternating signs in 
the (v,iS) space. This low-frequency divergence becomes 
more and more prominent when uj c increases steadily. In 
addition, one can distinguish two "cone" -like structures 
easily which are totally absent when uj c < 4.0. Next we 
pay attention to the Rel^ u part (see bottom panel in 
Fig El). The most striking features in the contour plots 
is the cross structure. Interestingly, the vertices can be 
split into two types based on the intensity and sign of the 
cross structures as well, (i) For oj c < 4.0, the cross struc- 
ture is notable and positive. Upon increasing the screen- 
ing frequency u> c , the intensity weakens and the extent 
shrinks, (ii) For oj c > 4.0, the cross structure exists too, 
but its sign is now inverted as indicated by the colors. In 
other words, the cross structures are concave, contrary 
to those in the uj c < 4.0 case. Furthermore, the vertex 
tends to be featureless with the increment of u> c . Thus, 
it is suggested that the frozen moment phase with large 
local magnetic moment will strengthen the two-particle 
vertex functions T, which is in accord with the previ- 
ous results obtained by Hafermann et al. [2l| They found 
that the vertex is essentially featureless on the Fermi liq- 
uid side, and evident in the frozen moment phase. To 
make a long story short, dramatic changes will occur in 
both the intensity and structure of Rex^ u and Rer^ u , 
accompanying with the first-order metal-insulator tran- 
sition and collapse of local magnetic moment, when the 
screening frequency uj c approaches the critical point. We 
also make detailed analysis on the other components or 
channels for x an d T, and the other reducible and irre- 
ducible vertex functions. [22| They exhibit analogous fea- 
tures when the dynamical screening effect is activated. 

(c) Dynamical screening effect on local two-particle 
vertex functions in SrVO^. So far all the discussions 
are limited to the simplest single-band Hubbard model, 
and the frequency-dependent interaction U is modeled 
by analytical Ohmic model. Q In order to confirm the 
above conclusions, we conduct further DFT + DMFT 
calculations for strongly correlated metal SrVC>3 to ex- 
amine the consequences of dynamical screening effect on 
the two-particle Green's functions and vertex functions. 
SrVC>3 is an ideal system to benchmark the new theories 
and tools describing the strongly correlated materials. 
It has been the subject of numerous experimental and 
theoretical studies. [3j In previous works, several DFT 
+ DMFT calculations were devoted to study the one- 
particle quantities of SrVC>3 under the dynamical screen- 
ing effect. @, 12, EH It is shown to renormalise the 



spectral weight near the Fermi level, to increase the ef- 
fective mass, and to suppress the t2 g quasiparticle band- 
width distinctly. Nevertheless, to our knowledge, the 
two-particle quantities of SrVC>3 have not been studied 
before neither experimentally nor theoretically. So our 
results can be regarded as a valuable prediction. 

The frequency-dependent interaction for realistic ma- 
terials can be evaluated at the random phase approxi- 
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FIG. 4. (Color online) Real part of the spin-up-up compo- 
nent of vertex functions for strongly correlated metal SrVC-3 
obtained in the DFT + DMFT calculations, (top) Local 
two-particle Green's functions Rex^'"- (bottom) Local ir- 
reducible vertex functions ReFj rr . The second fermionic fre- 
quency v' is fixed to be nT, and the bosonic frequency w 
is set to be zero. The static U is 4.0 eV, and the data for 
frequency-dependent U are taken from reference [f| directly. 
The DFT + DMFT calculations are done for /3 = 10 and 
/3 = 30, respectively. 



mation (RPA) level as usual. As for SrVOs, the U(cu) 
data are taken from reference [B| directly. The details 
for our DFT + DMFT calculations can be found in ref- 
erence [l4[. The main calculated results, real part of 
spin-up-up component of local two-particle Green's func- 
tion Re%^/ W , are shown in the top panel of FigQJ Com- 
paring the calculated results with and without dynam- 
ical interaction, the difference is rather significant, es- 
pecially for the v — > region. Moreover, the difference 
is amplified at lower temperature. When the dynam- 
ical screening effect is included, the absolute values of 
Re%^ " are strongly reduced. Similar or opposite trends 
are found for the other two-particle vertex functions, such 
as the local irreducible vertex functions Fi rr (see the bot- 
tom panel of Fig|4]). The vertices r irr can be evaluated 
straightforward through the famous Bethe-Salpeter equa- 



tion: Fir,. = /3(xo 1 ~X 1 ): where the polarization bubble 
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Xo is computed from the fully interacting one-particle 
Green's function. [22j 

Finally, it should be stressed that in many strongly 
correlated materials, the dynamical screening effect 
is conspicuous. [1, [H| Werner et al. 15] have demon- 
strated that the fine electronic structures of hole-doped 
BaFe2As2 can be fully understood only if the dynamical 
screening effect is taken into account for the first time. 
Recently, Park et al.^llfy have studied the magnetic ex- 
citation spectra in BaFe2As2 through a two-particle ap- 
proach within standard DFT + DMFT framework. Their 
calculated results, including the local irreducible vertex 
Tirr, dynamic magnetic susceptibility x(q, iu>), and dy- 
namical structure factor S(q,ui), roughly reproduced all 
the experimentally observed features in inelastic neutron 
scattering. 130J Since in their calculations just the static 
interaction was used, we believe that incorporating the 
dynamical screening effect into the DFT + DMFT cal- 
culations will contribute to improving their results. 

In summary, we study several typical one-particle and 
two-particle quantities in both the half-filled single-band 
Hubbard model and transition metal oxide SrV03 by 
applying the DMFT and DFT + DMFT methods with 
frequency-dependent interaction U. In the case of the 
Hubbard model, the dynamical screening effect is shown 
to drive first-order metal-insulator transition and mag- 
netic moment collapse by tunning the screening parame- 
ter uj c of the Ohmic model. The intensity and structure 
for both the local two-particle Green's functions Xaa'^ 
and full vertex functions T^," display dramatic changes 
when the critical w c is attained. Concerning the strongly 
correlated metal SrV03, the calculated two-particle ver- 
tex functions are also suppressed or enhanced by the dy- 
namical screening explicitly. Enormous changes for the 
q-dependent response functions are expected. Anyway, 
the two-particle quantities in either Hubbard model or re- 
alistic materials are intensively affected by the dynamical 
screening effect. Hence, in future studies with regards to 
moment-resolved response functions and diagrammatic 
extensions of DMFT, careful considerations for the dy- 
namical screening effect should be essential. 
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